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ABSTRACT 


This work investigated the suumertime northwestward propagating waves 
over the tropical western North Pacific and their relationship with tropical cyclone 
activities. The 15-season (May-October 1974-1988) Navy tropical global band data 
were separated into three groups according to the degree of organization in the 
multiple-set canonical correlation (MCC) modes that were computed from the 
surface v wind. Groups A (6 seasons) and B (5 seasons) showed orderly 
large-scale propagating patterns at the surface and 700 hPa with an average period 
of 8 days and a wavelength of about 2500-3000 km. The structure is robust 
whether the MCC analysis was performed on each of the seasons individually or 
on groups of years. A strong in-phase relationship between the low-level 
disturbance cyclonic cells and tropical cyclone centers was found. Composite 
analysese of all cases in these two groups whose MCC mode 1 amplitude exceeded 
1.5 m/s showed organized northwestward propagating divergence/convergence 
patterns at 200 hPa. Group C (4 seasons) shows the least organized structure, with 
a large difference between MCC modes obtained from the entire group and those 
from each individual season. The possibility that group C contains irregular 
tropical cyclone tracks was discussed. 
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I. INTRODUCTION 


The tropical easterly wave disturbance was first 
recognized by Piersig(1936) and Regular(1936), using surface 
data, they found these disturbances traveling from western 
Africa toward the Atlantic. These waves disturbances have a 
relatively well defined period of around 4 or 5 days with a 
horizontal wavelength of 2000-4000 Km in the lower 
troposphere. Dunn(1940) used surface pressure data to trace 
the movement of an isallobaric center across the western 
Atlantic ocean and Caribbean sea. During and after World War 
II ,using upper air data from the Caribbean area, 

Riehl(i945) and associates developed the theory of easterly 
wave model. The term " easterly waves " is used to described 
these synoptic disturbances, and their propagation 
characteristic are believed to be near those of Rossby 
waves. These waves, which have maximum intensity in the 
layer from 700 hpa to 500 hpa, slope eastward with height. 

Yanai et al. (1968) also found westward propagating 

waves, which have 4-5 days period with 6000-10000 Km wave 
length and eastward tilt vertical structure. Between the 
late 1960's and the early 1970's, most of the study of 
tropical tropospheric synoptic disturbances, using station 
radiosonde data and satellite cloud pictures, were made over 
the tropical western pacific. Over the western North 
Pacific, Wallace and Chang 1969, Chang et al. 1970 and Reed 
and Recker 1971 used spectrum and composite analyses to 
study the characteristics of the easterly wave disturbance. 
They found that these waves have a well-defined 4-5 day 
period and with a horizontal wavelength of 2000-4000 Km 
there is also a northeast-southwest horizontal tilt, which 
implies by a northward momentum flux. The waves contain a 
warm core vertical structure, which is equivalent barotropic 
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with very little tilt the over central part of the western 
Pacific, and have a slight westward tilt with height over 
the western part of the western Pacific. The influense of 
the vertical shear in the ambientflow, on the change in 
vertical tilt across the western Pacific bases was verified 
in a linear model by Holton (1971). 

The easterly waves are Rossby like waves. Over the 
eastern part of the western Pacific, the disturbance 
sometimes have a longer-wavelength. The vertical structure 
contains an eastward tilt with height below 300 hpa and 
westward with height above it. These waves are sometimes 
identified with mixed-Rossby-gravity waves. Using the 
internal variation of sea surface temperature Chang and 
Miller (1977) and Dunkerton(1993) and Chang and Zambresky 
(1994) , used suced seasonal variation to explain the change 

in the wave disturbance. 

The expanded observation platforms of the 1979 First 
GAPP Global Experiment (FGGE) facilitated further 
investigation into the behavior of tropical wave disturbance 
at locations where only few station data were available 
before. Nitta et al. (1985), Nitta and Takayaba (1985) and 
Tai and Ogura (1987) use these data to examine synoptic 
scale disturbance throughout the tropics. They identified 
various regions with active summer-time transient activities 
in the northern tropics by examining band-pass - filtered data 
of meridional wind(v) and outgoing longwave radiation (OLR). 
They concentrated on spectral characteristics and the 
vertical structure of the wave disturbance. Nitta and 
Takayabu (1985) and Tai and Ogura (1987), found west- 
northwestward propagating disturbances in the eastern 
tropical Pacific. They also observed cold anomalies below 
700 hpa and stronger warm anomalies at about 300 hpa. 
Recently, because of the increase in quality of numerical 


2 



model analyses and satellite data over the tropics, we are 
no longer limited to only station data. Reed et al. 

(1988a,1988b) studied the synoptic evolution of Africa wave 
disturbances appearing in the European Center for medium 
Range Weather Forecast (ECMWF) for the summer of 1985. They 
found two preferred tracks for Africa wave disturbance and 
also confirm much of the previous studies of the wave 
disturbances. They found over data sparse area the African 
wave disturbance features depicted by the ECMWF analysis are 
substantiated by independent satellite data. Liebman and 
Henton (1990) used the initialized data from ECMWF and 
outgoing longwave radiation (OLR) data during fall(September 
- December), 1980-1987 to study the zonal variation of the 
3- day equatorial synoptic scale waves. Lau and Lau (1990) 
also used the ECMWF data and OLR data for summer (June - 
August ), 1980-1987 to study the tropical waves disturbance 
over the western Pacific. They used single point correlation 
and rotated " extended empirical orthogonal function " 

(EEOF) analysis to show a northwestward phase propagation, 
with a wavelength around 2800 Km and period of 6 days. This 
result is within Liebman and Henton’s (1990) estimate. The 
summertime phase propagation pattern with the northeast- 
southwest elongated wave structure is more well defined than 
the fall pattern and it may be track back to 5'N,155-160‘E. 
Their results show that the mixed Rossby Gravity waves(MRG) 
will dominate the equatorial eastern Pacific and Rossby 
wave-like, west-northwest propagating waves dominate the 
western Pacific. The strong northwestward tilt with height 
is found in the equatorial region, and convection ( as 
indicated by OLR ) is ahead of low level troughs. The 
vertical tilt and convection-trough phase difference both 
decrease downstream, and disappear toward the south coast of 
China. 
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Takayaba and Nitta (1993), used the summer ECMWF data 
(1980-1989) along with satellite cloud data from Japanese 
GMS, to study the northwestward propagating wave pattern in 
the western North Pacific. They found the same westward 
propagating mixed Rossby-gravity wave (3-5 days) patterns 
which may dominate because to the mean vertical shear and 
sea-surface temperature. Their results in agreed with Chang 
and Miller’s (1977) finding,which was based on radiosoude 
data. They also suggest that these northwestward propagating 
wave disturbances may develop from the mixed Rossby Gravity 
waves in the equatorial central pacific. This downstream 
development appears plausible as both type of waves have the 
same 3-5 days period in their study. Nevertheless, Lau and 
Lau (1990) considered this link unlikely owing to the 
weakening of the central Pacific equatorial wave near the 
dateline. 

Chang et al.(1994) Who used multiple canonical 
correlation analysis(MCCA) with analyzed data from the Navy 
Operational Global Atmosphere Analysis and Prediction Model 
(NOGAPS) studied these northwestward propagating waves 
during the summer of 1989-1991. To focus in the synoptic 
scale disturbance they used high-pass filtered deta to 
remove features with a period greater than 12 days. The wave 
disturbance have patterns that are similar to those found by 
Lau and Lau (1990) but have a longer period of 8 days. They 
concluded that most of these waves are related to tropical 
cyclone activities, and postulated that nonlinear barotropic 
effects, which cause anticyclogenesis southeast of the main 
cyclone center with a southeast propagating of wave energy, 
favors a continuous development of the wave disturbance. 

These are Rossby-like waves where the beta effect plays a 
major role in the propagation and development. 

In addition to synoptic waves discussed above. 
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intraseasonal oscillations have been known to affect 
tropical weather significantly. Madden and Julian 
(1971,1972), using spectra and cross spectra with data from 
other stations, found two propagating waves , a westward 
propagating was found with a 5-6 days period and a global 
horizontal scale, ( Maddan and Julian 1972b ), the other is 
the 40-50 dsys oscillation which called Madden-Julian 
Oscillation(MJO). These 40-50 days waves propagate eastward 
with planetary scale Zonal wavenumber one. The dominant 
signal is in the zonal (u) wind component. There are many 
studies on the MJO and its’ role in tropical convection, 
especially over the Indian monsoon region(Krishnamurti and 
Subrahmanyan, 1982; Lorenc, 1984; Krishnamurti, 1985; Lau 
and Chan, 19 86) . 

Krishnamurti et al. ,(1985) clearly identifies a 
planetary-scale divergence wave that travel eastward around 
the globe throughout the FGGE year. The divergence or 
convergence centers usually show up over the western 
equatorial Indian ocean. These centers then propagate to the 
western Pacific. Hayashi and Sumi(i986) and Lau and 
Peng(1987), Used a numerical study to suggest that cumulus 
convection is an important factor in generating the 40-50 
day low frequency mode. The northwestward propagation of the 
divergent mode might be attributed to the supply of 
available water vapor and cumulus convection over the entire 
Asian Monsoon region. The theoretical studies have suggested 
that the MJO is likely moist Kelvin waves(Chang, 1977; Lau 
and Peng, 1987,-Chang and Lim,i988), and also may have a 
Rossby component that couple the convection with the 
planetary boundary layer(Wang, 1988). 

There is another intraseasonal signal that appear in 
the tropical wave disturbances that occurs in the period 
range between the synoptic wave and the MJO. Wallace and 
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Chang(1969) found that the meridional(V) winds at the 
tropical western Pacific show a 10-20 day peak. This quasi - 
biweekly periodicity was also found in various studies of 
the Asian summer monsoon region (Krishnamurti, 1973; 

Murakami and Frydrych, 1974; Zangvil, 1975; and Chan and 
Chen, 1993). These 10-20 day oscillations generally 
propagate westward and have a near barotropic vertical 
structure. Many of the previous studies associated this 
oscillation with the activity of the Indian summer monsoon, 
the studies of Krishnamurti and Bhalme (1976) and LI and 
Zhou (1992) show that the biweekly oscillation dominates the 
fluctuation of wind and surface pressure over the Indian 
monsoon region and documented its interaction with MJO. Harr 
(1993) studing the effect of variation of large scale 
circulation on tropical cyclone characteristics over the 
tropical western North Pacific. He identified six recurrent 
700 hpa circulation patterns that represent large-scale 
variabilities associated with the monsoon trough and 
subtropical ridge. There are many possible transitions among 
recurrent circulation patterns, but significant transitions 
occur over a very limited set of paths which are associated 
with interactions with tropical and midlatitude circulation 
systems. Furthermore their transitions can explain much of 
the observed intraseasonal variability in the occurrence and 
track types of tropical cyclones. 

The quasi-biweekly oscillations and the MJO both have 
strong signals over the northern summer monsoon region. This 
suggests that the intraseasonal oscillation may be the 
manifestation of a large-scale monsoon oscillation over the 
entire Indian Ocean-western Pacific region. (Krishnamurti 
and Subrahmanyan, 1982; Lau and Chen, 1986; Chan et al. , 
1988). The near 10 day period disturbance (Chang et al. , 
1994), nearly overlaps with the short end of the period 
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range of the quasi-biweekly oscillation. Most previous 
studies have relied on a time frequency filter, typically 
band-pass filters, to separate the MJO and quasi-biweekly 
systems. This may pre-determine the period window to be 
studied and may leave out structures that are outside of the 
predetermined window. The MCCA (Chen and Chang, 1993; Chen 
et al. ,1994) effectively isolate the synoptic scale, west- 
northwestward propagating disturbance in the tropical 
western Pacific from the NOGAPS analysis during summer of 
1989-1991, Although intraseasonal frequencies were removed 
by a high-pass filter that truncates the spectrum at a 
period of 12 days, narrow band pass filtering was not used. 
Because the technique of MCCA is mainly a spatial 
decomposition, is used to find the structure of the synoptic 
waves, which have a period close to 8 days during the three 
summers. This is longer than the previously assumed 3-6 day 
period. 

The purpose of this work is to study the properties of 
the tropical northwestward propagating wave disturbances 
over the western Pacific during 15-years(May-October, 1974 - 
1988), and to compare the results with tropical cyclone 
activities. The focus of this study will be on intraseasonal 
oscillations other than the MJO. So the MCCA will be applied 
to the meridional wind component. This will bring out the 
non MJO signal. Since the MJO is base on the zonal wind 
component. As we expected. The MCCA method successfully 
separates out tropical synoptic wave disturbances which we 
are interested in this study. 

In section 2 we describe the data and methodology used. 
Section 3 will show the leading MCC structure of v wind 
component. Section 4 will show one-point correlation and the 
relation between the synoptic wave disturbance and the 
tropical cyclone position. Section 5 discuss the results of 
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a composite study which show the contour and vector map of V 
vector from selected strong cases. Section 6 include a 
summary and the conclusion. 
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Figure l:Map of data area (showing core domain for MCCA and 
large domain for single point correlation and composites) 
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II. DATA AND METHODOLOGY 


A 15 year (May-October) wind data set is used, except 
for 1974, which only contain July-October. analyzed fields 
are from the operational Global Band Analysis (GBA) of the 
Fleet Numerical Meteorology and Oceanographys Center, 
Monterey, California, during 1974-1988. These data were 
produced four times daily by objective procedures on a 
mercator grid which extended from 60'N to 40‘S. The 
objective analysis scheme was designed to take advantage of 
all the reports from the operational data base. Although, 
data were available every six hours for the surface, 700, 
400, 250, and 200 hpa level. We only use twice daily data 
(00h and I2h) and surface, 700 hpa, 200 hpa (V) wind data. 
The domain of the study is shown in Fig l. A core region of 
110 -150 E, 10 S-25 N is embedded in the large domain of 
70 -180‘E, 20 *S-40’N. Because of major progress in computer 
technology and the availability of large amounts of data, 
researchers have used different kind of statistical methods 
in the meteorological data fields. One of these is the 
Principal Component Analysis (PCA) which was used by 
Pearson(1901) in a biological context when he recast linear 
regression analysis into a new form so as to avoid the weak 
asymmetrical relationship between "dependent" and 
"independent" variables. Lorenz (1956) used the technique in 
the field of statistical-dynamical approaches to weather 
prediction. The term PCA, sometime called "empirical 
orthogonal function analysis" (EOFA) concentrates the 
internal variation within one data field to a set of 
leading components. These are used to transform the original 
set of variables to a new set of uncorrelated variables 
called principle component(PC). These are used to find a 
linear combination of the original data variable such that 
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the principal component, which is a function of time over a 
set of grid point(17x15 in our study) and is also referred 
to as an "amplitude", has the maximum variance under the 
condition that the weighting function is orthogonal to those 
of all preceding PC’s. Subsequent PC’s can be found under 
the same maximum variance,unitary weighting function, and 
orthogonality requirements. If a data pattern is 
proportional to its own weighting function, all other PC’s 
will vanish due to the spatial orthogonality. These 
weighting functions are often called EOF’s in meteorological 
study, also we can reduce the dimensionality of the data in 
order to simplify later analysis. Chen et al., (1993) have 

shown that direct application of PCA to data from a well- 
chosen network is particularly efficient. 

PCA gives us the internal variation within one data 
field to the leading component, but does not tell us the 
information of relationship between different fields. The 
relationships between two data fields in term of correlation 
coefficient can be determined by the canonical correlation 
analysis (CCA). The CCA is basically a search for optimal 
correlation between two variable fields by linearly 
combining the elements of each fields. These combinations 
are called canonical component(CC). The CCA method pairs 
patterns (one for each field) so as to identify those 
parts of the patterns that are significantly correlated and 
their respective variance contributions. Conventional CCA 
can be used for investigating the correlation between two 
fields, but it can only show maximum correlations without 
regard to the variance significance. So we always perform 
the PCA before the CCA does. It is hoped that with PCA, the 
leading component may represent different physical 
processes, or at least, component truncation can help to 
delete noise and errors which may be associated with small 
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variabilities. 

The multiple-set canonical component analysis(MCCA) 
method is a generalization of the traditional two variable 
CCA to multiple variables(e.g. Steel, 1955; Horst, 1965; 
Kettenring, 1971; Gifi, 1990). This method tries to find the 
maximum correlation among more than two data fields. Through 
a diaganolization of the product or the square product of 
the correlation matrices, Chen et al.,(i994) not only 
simplify the MCCA procedures to be used in studies with 
large dimensioned data, but also can select the correlation 
between the two fields we need. The results are accepted 
when the largest residual (cross-component) correlation is 
negligible compare to the average correlation. The 
application of the MCCA method follows that of Chen and 
Chang(1993,1994), where the different data fields are 
constructed from an overlaping-sliding time window that 
contains consecutive fields of the same variable. In our 
study, we use MCCA on v s£c , 7 00 and 200 hpa within the core 
region to find the basic disturbance structure mode, the 
basic structure mode is correlated with other variables to 
find the entire structure of the disturbance. 

The data on 255 grid points within the core domain are 
first pre-processed with an EOF analysis, that will truncate 
the data to include the first 18 principal components, which 
retain 79% of the total variance in the core domain. A high 
pass filter, 20 days, is applied to each of the components 
by removing all waves whose period longer than 20 days, 
which may masks significant synoptic time scale events. This 
will leave the higher frequency, synoptic-scale events. The 
data set, which contains 12 consecutive twice-daily data, 
are formed as the window (00h(time to) to 132h(time 
til)slides forward at 12-hour intervals. The squared product 
of correlation matrices between fields of adjacent times are 
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optimized. The reliability of this application procedure of 
the MCCA method for the field data has been tested and found 
satisfactory by Chen and Chang(1994). 

The MCCA results in 12 MCC modes that are ranked 
according to the average correlation through any 132 h(5.5 
days) period. The first mode always contains the largest 
correlation coefficience, which will show the best defined 
pattern. The comparisons between the wave patterns in 
different time frames will tell the characteristics of 
transient phenomena occurring in the area covered in this 
study. 

At any given time MCC mode contains both horizontal 
spatial and time structure, which includes the 12 
consecutive 12-hour time frame from day 0 to day 5.5. 
Therefore, we also produce a time-integrated amplitude 
coefficient. This coefficient is the average amplitude of 
the 12 consecutive 12-hour structure and is plotted in the 
time domain according to the first time frame,(tO). 

To study the U,V structure over large domain and at 
other level, we use one point correlation. The v sfc , v 700 , v 200 
wind component are determined by the single point 
correlation between these field and the time-integrated 
amplitude coefficient of the MCC mode one. A sequence of 
correlation plots based on different time lags may be 
produced. The correlation plots will give us some idea of 
how well the MCC mode is done in the mean removed data. We 
also use the correlation plots to find the possible 
relationship between the propagating wave and tropical 
cyclone activities. This is done by identifing the time of 
the maximum point in the time series and then correlate it 
with the tropical cyclone center positions. 

Another approach to reveal the three dimensional and 
time pattens of the leading MCC modes are to composite the 
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strongest cases chosen on the time-integrated amplitude 
coefficient time series. The composite was with respect to 
the time lag from the large amplitude time points. Because 
the cases selected are based on a particular MCC mode, it is 
likely that each case is indeed similar and that the 
composite does provide a useful representation. 
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III. LEADING MCC STRUCTURE 


All MCC modes are determined from the May-October v Bfc 
fields. Several calculations are performed based on 
different periods of the data. First, the entire 15-year 
summer data are used to find a set of MCC modes, which are 
referred to as the 15Y MCC modes. The MCCA is then 
performed on each of the 15 individual years, resulting in 
15 sets of MCCs, one for each year. The individual year's 
leading MCC mode are compared to the 15Y's, and are found to 
resemble the 15Y leading mode in various degrees. We then 
group the individual years into three groups, with Group A 
containing years whose MCC mode 1 resembles that of 15Y, 
Group B contains years whose MCC mode 1 resembles 15Y to a 
slightly lesser degree, and Group c containing years whose 
MCC mode 1 dose not resemble that of 15Y. Finally, the MCCA 
is performed on each of the three groups, resulting in yet 
another three sets of MCCs. 

In most cases, MCC modes 1-2 form a pair of propagating 
modes with similar horizontal structure, fractional 
variance, averaged correlation, and approximately one 
quarter cycle phase difference. Higher order modes are 
normally associated with small variance and smaller scale, 
less coherent structure. Thus, our study focuses on mode 1 
only. 

A. MCC MODES DETERMINED FROM 15-YEAR DATA 


The first two MCC modes for the 15Y are shown in Fig. 

2. The 12-hour frames show a synoptic-scale, 
west-northwestward propagating pattern with a wavelength of 
about 2500-3000 km. The propagation track starts from the 
tropical western Pacific, across the Philippines and reaches 
southeastern coast of China. Any two time frames that are 
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96 h apart demonstrate a near 180* out-of-phase relationship 
in the main features. For example, time frame T2 shows a 
positive cell center over the northern Philippines and a 
negative cell center at 15*N, 135*E; while time frame T10 
shows a negative center over the northern Philippines and a 
positive cell center at 12*N, 133‘E. Thus, a period of 
approximately 8 days may be inferred. These horizontal 
structures and periodicity are very close to those estimated 
by Chang et al. (1995) based on NOGAPS data of 1989-1991 

using the same MCCA technique. The structures also resemble 
those found by Lau and Lau (1990) using simple lag 
correlations on the 1981-89 ECMWF data, although they 
estimated a period range of 6-8 days. 

The averaged correlation between each consecutive 
12-hour frame is about 0.87 for both modes 1 and 2, which is 
slightly less than the 1989-1991 NOGAPS results of 0.94. 

The total fractional variance for the two modes is about 
25%. Because of different filters used, a direct comparison 
of the present results and the three-year NOGAPS study can 
not be made, but it is possible to state that this variance 
is slightly less than that obtained by Chang et al. (1995). 

These slightly reduced correlations and variances are not 
significant, considering the fact that the GBA data set does 
not use model forecasts as a first guess. The present 
variance estimate is significantly larger than Lau and Lau's 
(1990) results of about 9% for the northwestward propagating 
wave. This is likely due to the bogus typhoon data used in 
the GBA which is not used in the ECMWF analysis. In 
addition, Chen and Chang (1993) and Chen et al. (1994) found 
that the multiple-time correlations over a synoptic period 
can isolate a significantly larger signal than simple lag 
correlations. The latter represents only a partial view of 
the time variability based on two time points. The signals 
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isolated depend on how much these two time points can 
represent the entire variability. On the other hand, the 
former method seeks for the variability based on a more 
complete sampling of several time points throughout the 
synoptic period, and may give rise to a more complete 
representation of the coherent motion system. 

B. MCC MODES DETERMINED FROM INDIVIDUAL YEARS 

A main purpose of the present study is to see how 
representative the 8-day northwestward propagating wave is 
for each individual year in the GBA data set We therefore 
conducted an MCCA for each of the 15 years. In all cases 
modes after the first two do not contain consistent-varying, 
large scale structures throughout the 5.5 day period of the 
MCCA, so again only the leading pair of MCC modes 1-2 are 
considered interesting. The individual MCC mode 1 
structure, derived from individual year from 1974 to 1988, 
are included in Appendix A (Fig. Al). In general, the 
structures may be classified into three broad categories 
based on their resemblance or difference from the 15Y mode 
l. These categories are termed groups A, B, and C, 
respectively. Group A contains the six years that best 
resemble the 15Y. Group B contains five years that partially 
resemble 15Y. Finally, group C includes four years that 
appear to be significantly different from 15Y. The years 
within each group are listed below: 

Group A:1975, 1977, 1979, 1985, 1986, 1987 
Group B:1974, 1980, 1981, 1982, 1983 
Group C:1976, 1978, 1984, 1988. 

In this categorization Group B years usually show the 
basic cell structure of the 15Y pattern, but are considered 
to have only partial resemblance because either the 
structure in a number of time frames shows smaller - scale or 
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other anomaly features, appreciably depart from the straight 
northwestward track, or the whole pattern appears stationary 
over part of the period. 

The total fractional variance and the averaged 
correlation for MCC modes 1 and 2 for the individual years 
are listed in Table I. The fractional variance, in general, 
ranges between 24% and 40%, mostly above the total variance 
determined from the 15Y MCC modes 1-2. The single exception 
is 1978, a group C year, with a very low variance of 8%. We 
have examined other modes of this year, none of which show 
patterns that resemble the leading 15Y MCC modes. The 
variance of group A years, on the average, is slightly 
higher than the other groups, and that of group C years, 
even after removing 1978, is the lowest, although the 
differences are not large. The correlation for individual 
years ranges between 0.75 and 0.93, as compared to the 15Y 
value of 0.87. There is no systematic change among the 
three groups. 


20 




YEAR 

GROUP 

TOTAL 

VARIANCE 

CORRELATION 

15 


25% 

0.87 

1974 

B 

37% 

0.91 

1975 

A 

27% 

0.89 

1976 

C 

31% 

0.93 

1977 

A 

31% 

0.90 

1978 

C 

8 % 

0.90 

1979 

A 

40% 

0.93 

1980 

B 

33% 

0.93 

1981 

B 

32% | 

0.88 

1982 

B 

28% 

0.80 

1983 

B 

29% 

0.85 

1984 

C 

26% 

0.88 

1985 

A 

39% 

0.92 

1986 

A 

38% 

0.89 

1987 

A 

34% 

0.85 

1988 

C 

24% 

0.75 


Table I:Variance and correlation for 15Y 
and individual years MCC mode#l and #2 

Since individual years in group have their leading MCC 
mode patterns appreciably different from the 15Y, it will be 
interesting to see how an MCCA for the entire group 
(combining all years within the group) will compare with the 
15Y patterns. Therefore, we carried out the MCCA for each 
of the individual groups. The resultant MCC modes 1-2 
patterns are shown in Figs. 3, 4, 5, respectively, for 
groups A, B, and C. The total fractional variance and the 
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average correlation for these group modes are shown in Table 
2. As expected, the group A mode 1 pattern (Fig. 3a) shows 
the northwestward propagating wave pattern with similar 
spatial and time scales as that of 15Y. In this case mode 2 
of group A (Fig. 3b) is nearly in phase with the 15Y mode 1 
(Fig. 2a), and vice versa. The total modes 1-2 fractional 
variance is 32% (Table II), which is considerably higher 
than the 15Y value of 25%. The average correlation of 0.88 
is very close to the 15Y's 0.87. 


YEAR 

GROUP 

TOTAL 

VARIANCE 

CORRELATION 

A 


32% 

0.88 

B 


28% 

0.86 

C 


13% 

0.83 

CC 


23% 

0.83 


Table II:Same as Table I except for Group 


The group B modes 1-2 patterns (Fig. 4) also resemble 
the 15Y patterns, almost as close as that of group A, with 
its mode 1 corresponding to the 15Y mode 2 and vice versa. 
The resemblance is an improvement over those of the 
individual B years. Therefore, the five years, when 
combined, give a more clearly defined large-scale wave 
signal than individual years. The total fractional variance 
is 28%, and the average correlation is 0.86, which is again 
very close to the 15Y's value. The group C modes 1-2 are 
shown in Fig. 5. Different from each of its four individual 
years. Fig. 5 reveals a pattern that contain some aspects of 
the northwestward propagating wave pattern shown in the 15Y 
modes. Fig. 5b shows that main features of group C's mode 2 
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have a small phase difference (somewhat less than 90’) with 
the 15Y mode l (Fig. 2a). In the meantime, the main 
features of group C's mode 1 (Fig. 5a), particularly in the 
vicinity of the Philippines, are roughly half a cycle 
out-of-phase with the 15Y mode 2 (Fig. 2b). Thus, group C 
modes 1-2 roughly describe a propagating wave that is about 
180‘ out-of-phase with the 15Y modes 1-2. This partial 
resemblance of the 15Y modes by the group C result is a 
significant difference from each of the four individual 
years. It indicates that the common structure among these 
four "abnormal" (relative to the 15Y structure) years is 
closer to the 15Y modes than that is revealed by any 
individual year's MCC modes. The total mode 1-2 fractional 
variance of group C is 13%, considerably smaller than the 
15Y, apparently because of the very low variance in 1978 
(Table l). The average correlation is 0.83, which is only 
marginally less than 15Y and groups A and B. 

In order to see whether the low-variance 1978 has any 
special impact on the group C results, we conducted the MCCA 
on the three other years of group C. This shall be termed 
group CC, and the resultant MCC modes 1-2 are shown in Fig. 
6 . The patterns are much the same as those from group C. 
Table II also shows that the total mode 1-2 fractional 
variance for group CC is 23%, considerably higher than group 
C (when 1978 is included in the group), but smaller than 
group A's 32% and group B's 28%. This systematic decrease 
in the fractional variance from groups A to C suggests that 
the northwestward propagating wave pattern is the domainant 
organized motion system in this region during the northern 
summer . When this motion system is more active, the v s£c 
field is more organized, which leads to larger fractional 
variances of MCC modes 1-2. The fact that groups B and C's 
modes reveal a structure closer to the 15Y modes than the 
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individual years within their respective groups is an 
additional indication that the northwestward propagating 
wave pattern is a common structure among most years. 

C. TIME VARIATIONS OF THE MCC MODE AMPLITUDES 

The time variations of the MCC mode amplitudes are 
shown both as a composite time series and a power spectrum. 
The reason of a composite time series is that, at any given 
time, a mode contains both horizontal and time structures, 
with the latter encompassing 12 consecutive twelve-hourly 
frames from day 0 to day 5.5. Thus, the time variation of 
each mode should be described by considering the fact that 
each mode is a space-time volume function that covers not 
just the horizontal spatial domain but also a time domain of 
5.5 days. In order to take into account this time domain, 
the coefficient of each mode is constructed by a time- 
integrated amplitude coefficient. This coefficient is the 
average amplitude of the 12 consecutive twelve-hour 
structures, and is plotted in the time domain according to 
the time of the first (day 0) time frame. 

we have conducted the MCCA on the 15-year data set as 
one group (15Y), in three different groups (A, B, C), and 
for each individual years, For each year three 
time -integrated amplitude series of mode 1 can be 
constructed. Fig. 7 shows the three time series for each of 
the six group A years, and Figs. 8-9 show those for the 
groups B and C years, respectively. For each year the top 
panel is the time series from the 15Y mode 1, the middle 
panel is from the respective group's mode 1, and the lower 
panel is from mode 1 based on the single-year data. The 
power spectra are shown to the right of the time series. In 
general, the six member of group A show a high correlation 
among the three time series, particularly between the 15Y 
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and the group A series, which are mostly in-phase throughout 
the season. The single-year series are also correlated to 
the respective 15Y and group A series, an occasional a phase 
shift of around one quarter cycle is simply due to 
identification of the TO time frame in applying the MCCA 
procedure to different lengths of data. As long as the mode 
pattern also has a one-quarter cycle shift (approximately 
two days), the results coincide with each other. In some 
years (e.g., 1985), the single-year phase shift during the 
early part of the season (May-July) is larger than during 
the late part of the season (August-October). As will be 
shown in the next section, the northwestward propagating 
wave patterns can often be associated with tropical cyclone 
activities. These activities are normally more active 
during the latter period. Therefore, the better 
correspondence during late summer between the single-year 
MCC mode l series and that of the group A and 15Y may be due 
to the increased presence of tropical cyclones. 

The good correspondence among all three time series for 
each year is also obvious in the group B years (Fig. 8) 

Since the group C modes 1-2 describe a propagating wave 
pattern that is roughly 180' out-of-phase with 15Y, in Fig. 

9 the sign of the group C series has been reversed in order 
to help the comparison. For group C years, the 15Y and 
group C time series are in good agreement, but their 
correspondence with the single-year series is less. This is 
consistent with the fact that the single-year MCC modes 
least resemble the 15Y modes. 

Figs. 7-9 show that the power spectra for the different 
MCC modes and different years often have spectral peaks in 
the 6-12 day range. The most frequent peak occurs around 8 
days. In order to identify the dominate spectral peaks, we 
composited the spectra shown in Figs. 7-9 for the entire 15 
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years, and for each of the three groups. The 15-year 
composite is based on the 15Y modes 1-2 spectra and is shown 
in Fig. 10a. Each group composite is based on the group's 
modes 1-2 (not the 15Y modes), and is shown in Figs. lOb-d 
for groups A, B and C, respectively. The 15Y composite 
shows a broad, but well-defined, power peak in the 6-12 day 
period band. The group A composite shows a sharply defined 
peak at around 8 days, the same period that was observed in 
the 15Y MCC modes 1-2 patterns (Figs. 2-3), and is the 
period that was reported by Chang et al (1995) using the 
northern summer 1989-1991 NOGAPS data. The group B 
composite is different from group A. It shows a double peak 
at 6 and ll days, respectively, and the 8-day period appears 
as a local power minimum within a broad 6-12 day maximum. 

The group C composite contains large difference from group 
A. There is a broad spectral maximum that covers a range 
from 6-days to longer periods. The drop in power beyond 
20-days in Fig. lOd is likely due to the use of a high-pass 
filter. Thus, group C, which has the least organized signal 
of the northwestward propagating waves in its MCC modes, 
also has the smallest percentage of its power in the 6-12 
day range. 
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Figure 2a:Weighting function for 15Y MCC mode#1(equivalent 
to MCC mode#l structure) of surface meridional wind(v) for 
12 consecutive 12-hour frames from 00 h to 132 h(5.5 day). 
Contour interval is 0.06 and dashed lined correspond to 
northerly winds when MCC mode#l amplidtude is positive 
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Figure 3b: Same as Figure 2 except for Group A MCC mode#2 
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Figure 7b: Same as Figure 7a except for year 1977 
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Figure 7d: Same as Figure 7a except for the year 1985 
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Figure 8b: Same as Figure 7 except for the year 1980 
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Figure 10a: Composite power spectrum for 15Y MCC mode#l and 
#2 combined. 
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Figure 10b: Same as Figure 10a except for Group A 
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Figure 10c: Same as Figure 10a except for Group B 
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Figure lOd: Same as Figure 10a except for Group 
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IV. ONE-POINT CORRELATIONS 


The spatial patterns of this mode over the larger 
domain (70’E-180*, 20'S-40’N) in the meridional (v) and 
zonal (u) wind components at three levels, sfc, 700 hPa, and 
200 hPa, are determined by the single-point correlation 
between these fields and the time-integrated amplitude 
coefficient of mode 1. Fig. 11 shows the correlation 
patterns of v sfc with the 15Y mode 1, at lags 0 and 4 days, 
respectively. Here all 15 years of data are used. Since 
the MCC modes are calculated based on v sfc , the correlation 
pattern at zero lag closely resembles the structure of mode 
1 shown in Fig. 2. 

Chang et al. (1995) demonstrated that the northwestward 

propagating wave patterns are correlated with occurrence of 
tropical cyclones for the 1989-1991 summers of NOGAPS data. 
Adopting their method, the possible relationship between the 
propagating waves and tropical cyclone activities may be 
examined by comparing the time series of the time integrated 
wave amplitudes and typhoon occurrences. From the top 
panels of Figs. 7-9, the times of large-amplitude wave 
activities for the 15Y mode l are determined to compare with 
tropical cyclone activities. This is done by identifying 
the time of the maximum and minimum points for several 
phases of the disturbances when the time-integrated 
disturbance amplitude exceeds 1.5 m/s and then correlate it 
with the location of tropical cyclone centers as reported by 
the Joint Typhoon Warning Center (JTWC), Guam. These 
locations of tropical cyclone centers at maximum amplitude 
times of mode 1 are plotted against the v sfc -mode 1 
correlation map at phases TO and T4 days (Fig. li). For 
each phase the majority of the cyclone centers are 
concentrated in a region to the west of positive cells. 
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where the main positive dv/dx wave axis as defined by the 
correlation pattern is located. The clustering of tropical 
cyclone centers near this axis indicates a strong 
relationship between the main wave pattern in the core 
domain and the circulations of the tropical cyclones. This 
relationship is clearest in the western Pacific east of the 
Philippines, but it also holds in the northwest corner of 
the core domain near the southern China coast and the 
northern South China Sea. 

Since the tropical Global Band data contain bogussed 
tropical cyclone vortices at the surface, it is conceivable 
that the good correspondence between the tropical cyclone 
locations and the wave pattern is due to the bogussed data. 
However, the wave pattern describes a much larger spatial 
scale than the bogussed vortices. The bogussed tropical 
cyclones are not applied to other levels, and the objective 
analysis scheme is two-dimensional only, so the correlation 
map for v 700 provides information that is independent of the 
bogus data. Fig. 12 shows the correlation map of v 700 with 
the 15Y mode 1 phases TO and T4, respectively, for the 
entire 15 years. It is readily seen that the northwestward 
wave pattern is approximately reproduced in this Figure, and 
is nearly in-phase with the v sfc maps. The maximum 
correlation coefficient is around 0.3-0.4, which is only 
slightly below the v sfc values. Therefore, the wave pattern 
is contained in the v 700 data and not just an artifact of the 
bogussing procedure. 

The v 2 oo ' 15Y mode 1 correlation patterns for TO and T4 
are shown in Fig. 13. The maximum correlation values of 
around 0.1 are very small. Nevertheless, the correlation 
pattern appears to be organized and consistent with the 
low-level patterns. The TO phase (Fig. 13a) appears to be 
one quarter cycle out of phase with the v s£c (Fig. lla) , with 
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the cells tilted forward (northwestward). The Tl phase 
(Fig. 13b) appears to be half a cycle out of phase with the 
v sfc (Fig. lib). The weak correlation and somewhat different 
phase differences from v sfc makes the v 200 correlation 
patterns only marginally useful, indicating only a 
possibility that there is an out-of-phase structure 
contained in the v 200 data. 

The 15Y correlation patterns and the tropical cyclone 
location relationships suggest that the 8-day northwestward 
wave patterns, during large amplitude fluctuations, appear 
to be a manifestation of the tropical cyclone activities. 

It would be interesting to see how this organization of the 
wave pattern may appear in the different groups that were 
categorized according to the single-year MCC modes. In the 
ensuing Figures, we compare the tropical cyclone locations 
with the timing of the large amplitude group MCC mode 1, and 
plot them on the v sfc correlation maps for each of the three 
groups. For each group, only the data within the years 
belonging to that group are used. 

Fig. 14 shows the v sfc - group A mode 1 correlation 
patterns, with the tropical cyclone locations plotted in the 
phase TO and T4 maps, respectively. Here the tropical 
cyclone locations are also lined up in the positive dv/dx 
axis in the western Pacific, although for the cell near the 
northwestern corner of the core domain, which partially 
covers the southern China coast, the distribution of the 
tropical cyclones is more scattered within the positive 
cell. The maximum correlation at v sfc is between 0.4-0.5, 
which is slightly higher than the 15Y results shown in Fig. 
ll. The group A v 700 correlation patterns are shown in Fig. 
15. There is again a nearly in-phase relationship with the 
v s£c , and the maximum correlations of between 0.4-0.5 are 
again slightly higher than the 15Y results. The group A v 200 
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correlation map (Fig. 16) also shows a weak but organized 
pattern, in this case both TO and T4 are half a cycle 
out-of-phase with v sfc . The maximum correlation values are 
0.1-0.2, again slightly higher than those shown for 15Y. As 
group A includes the six years with the most organized wave 
pattern, it is not surprising that the correlation values 
are slightly higher and that the v 200 patterns are more 
consistent. 

Fig. 17 shows the v sfc - group B mode 1 correlation 
patterns, with the tropical cyclone locations plotted in the 
phase TO and T4 maps, respectively. The results are similar 
to those of group A, with maximum correlation around 
0.4-0.5, and tropical cyclone locations concentrated in the 
positive dv/dx axis in the western Pacific. Near the 
northwestern corner of the core domain the tropical cyclone 
centers are again more scattered. The v 700 correlation 
patterns (Fig. 18) are again in-phase with v sfc , with maximum 
correlation values around 0.3-0.4, slightly less than group 
A and comparable to 15Y. The v 200 correlation maps (Fig. 19) 
still show a weak but organized pattern with maximum values 
between 0.1-0.2. This pattern is roughly out-of-phase with 
low levels in the vicinity of the maximum correlation. 
Therefore, the five years in group B show the same organized 
wave patterns as group A. 

The v sfc correlation patterns for group C are shown in 
Fig. 20. The patterns are not as organized as the other two 
groups and resemble the 15Y patterns the least. At phase TO 
(Fig. 20a), the maximum correlation is only 0.2, 
significantly lower than other groups. The distribution of 
tropical cyclone centers are also more scattered, although 
there is still some indication of preferred locations near 
the positive dv/dx axis. At phase T4, the maximum 
correlation reaches 0.4, which is comparable to 15Y, but the 
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tropical cyclone centers deviate more from the positive 
dv/dx axis, and are mostly scattered within the positive 
cell in the core domain. As found with other groups, the v 700 
correlation patterns (Fig. 21) show a nearly in-phase 
relationship with the major cells of v sfc , but the maximum 
correlation drops to between 0.2-0.3. The v 200 correlation 
patterns, shown in Fig. 22, completely lack organized wave 
cell patterns with randomly low correlations. Therefore, 
group c years show the least organized wave patterns. 

We have previously noted that within group C, the 1978 
single-year MCC modes 1-2 contain only 8% of the total 
variance, significantly below the 24-40% range of other 
years (Table l). In order to find out whether the poor 
organization of the wave pattern in group C is mainly due to 
this anomalous year, the correlation maps are plotted for 
group CC, which excludes 1978 from group C. Figs. 23-25 
show the group CC correlation maps for v sfc , v 700 , and v 200 , 
respectively. In general, the organization (or lack of) of 
the patterns are similar to the group C results and are 
inferior to the other groups. The only significant 
difference occurs in v 200 (Fig. 25) . Here the northwestward 
propagating wave pattern is still missing, but along and 
north of 25'N (the northern boundary of the core domain) 
there appears to be a wave pattern with alternating signs. 
However, this looks distinctively unrelated to the low-level 
patterns and the correlations are very weak (0.1) . 

Therefore, we can conclude that the four group C years are 
anomalous years, which do not show the clearly-organized 
wave patterns, and corresponding tropical cyclone 
activities, that appear in other years of the 15-year data 
set, 
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lla: Single point correlation patterns of vsfc with 
MCC mode#l, at lag 0 days and with the tropical cyclone 
position ploted. Contour interval is 0.1 and negative values 
are dashed 
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V. COMPOSITE ANALYSIS 


Another approach to reveal the three-dimensional and 
time patterns of the leading MCC modes is to composite the 
motion fields according to the phases of oscillations in the 
MCC mode time series. This is carried out by identifying 
oscillations in the time-integrated amplitude time series 
whose signals are considered strong. For each case selected 
the timing of maximum amplitude is assigned to be time TO. 
The time variation of all wind fields are then composited 
with respect to the time lag, at 12 h intervals, from TO, 
for all selected cases. The time-mean background flow is 
removed from the composite charts in order to make the 
anomaly patterns standing out. The selection of strong 
signal cases are identical to that done in comparing 
tropical cyclone locations with the v sfc correlated pattern 
carried out in the previous section. Here, the timing of 
each large-amplitude occurrence for the MCC mode 1 is now a 
case for TO. This procedure gives 90 TO cases for the 15Y 
mode. (It should be noted that not all 90 cases have 
corresponding tropical cyclone locations in Fig. 11 , some 
cases have none and some have more than one.) As will be 
shown in Table III, these 15Y cases include 37 cases from 
group A years, 33 from group B, and 20 from group C. 



15Y 

GROUP 

COINCIDING 

A 

37 

42 

32 

B 

33 

32 

30 

C 

20 

19 

0 


Table III.Composite case in 15Y, Group A,B and C 

Fig. 26 shows the v sfc composite. As expected, the main 

TO and T4 patterns resemble closely the 15Y v sfc - mode 1 TO 
and T4 correlation patterns (Fig. li). During the first 24 h 
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(T0-T2) there is an indication of an northwestward extension 
of the wave pattern into the Asian continent, which is not 
found in the correlation pattern. The northwest extension 
of the wave pattern disappears as the time lag increases 
beyond T2, but the main wave pattern inside and near the 
core domain stays organized and propagates northwestward 
throughout the 5.5 days (TO-Til) of the composite. Fig. 27 
is a vector plot of the composite v sfc and u sfc . The wave 
pattern is dominated by a cyclone-anticyclone couplet that 
propagates through the core domain. The 15Y v 700 composite 
and corresponding wind vector plots are shown in Figs. 28 
and 29, respectively. The patterns are consistent with the 
correlation pattern (Fig. 12) and are approximately in-phase 
with the surface plots. In addition to the main wave 
pattern, during T8-T11 there appears to be a cyclonic vortex 
in the midlatitude western Pacific east of Japan (Fig. 29). 
Although there is some indication of a similar signal at the 
surface, this system does not appear to be consistently 
present or move through the entire series, and there is no 
corresponding correlation in Fig. 12 to suggest any 
significance. 

The 15Y 200 hPa v and wind vector composites are shown 
in Figs. 30 and 31. The v 200 composite shows a noisier 
picture than the low levels, but the propagating wave 
pattern inside the core domain can readily be seen. The 
structure shows more northeast-southwest elongation than 
circular cell features. In general, the propagating 
patterns are more conspicuous than those indicated by the 
15Y v 20C correlation (Fig. 13) . The wind vector composite 
(Fig. 31) appears noisy, but significant large-scale 
divergent and convergent patterns are noticeable. In Fig. 

31 the main areas of divergence and convergence in and near 
the core domain are marked. We can see a general pattern of 
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northwestward propagation of these divergence/convergence 
areas and also the propagation of anticyclonic and cyclonic 
features. The divergence pattern suggests that the wave 
mode is associated with strong convection, which is 
consistent with the observed association with tropical 
cyclone activities. 

The composite analysis is also carried out for each of 
the three groups. For group A, 42 cases are selected. 

Table 3 shows that 32 of these cases are also selected in 
the 15Y analysis. The surface v and wind vector composites 
(Figs. 32) are very similar to the corresponding 15Y 
composites, showing clearly organized northwestward 
propagating wave patterns. In fact, the group A patterns 
are more clearly defined than the 15Y, which is not 
surprising as group A contains the six most organized years. 
In addition, the northwest extension of the wave pattern 
over the Asian continent noticed during T0-T2 at the 15Y v sfc 
composite (Fig. 26) is also detectable in group A. In fact, 
the signal of this extension shows up slightly at T0-T1 then 
disappears, to reappear after T6 and becomes quite visible 
during T9-Til. Thus, the group A composite suggests the 
possibility that the wave pattern can continue to propagate 
beyond the core domain into the region of Tibetan. This 
downstream propagation is also visible in the wind vector 
composite (Fig. 33). The 700 hPa group A composites (Figs. 
34-35) also clearly depicts the wave pattern and resemble 
the 15Y pattern. The 200 hPa group A composites are shown 
m Figs. 36-37. The wind vectors indicate a similar 
northwestward propagating divergence/convergence pattern to 
that found in the 15Y group. The strong wind speeds in the 
Group A pattern reconfirms that the group A years are the 
strongest, most organized years for this wave. 

The group B composites contain 32 cases, of which 30 
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are included in the 15Y composites. In general, group B 
composites at the surface (Figs. 38-39) and 700 hPa (Figs. 
40-41) show clearly defined wave patterns that are 
comparable to the 15Y results. At 200 hPa (Figs. 42-43), the 
patterns also resemble the 15Y pattern with slightly stonger 
wind speeds. 

The group C composites contain 19 cases, none of which 
overlaps with the 15Y cases. These are the four years with 
least organization, and, as is with the correlation results, 
the wave pattern is only detectable at the lower levels. 

The surface and 700 hPa composites are shown in Figs. 44-45 
and 46-47, respectively. There are indications of the wave 
pattern but the Figures are also filled with other noisy 
features. The 200 hPa composites are shown in Figs. 48-49. 
There are a number of rotational and divergent disturbance 
features in each of the time panels, but there is no 
discernible time-continuous, organized structure in the core 
domain that can be related to the northwestward propagating 
wave pattern. 
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Figure 28: Same as Figure 26 except for v700 
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Figure 34: Same as Figure 26 except for Group A v700 
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Figure 36: Same as Figure 26 except for Group A v200 
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Figure 45: Same as Figure 27 except for Group C vsfc 
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VI. SUMMARY AND CONCLUSIONS 


This research used the MCCA method to study the 
tropical summertime synoptic-scale disturbances over the 
northwestern Pacific that were contained in the 1974-1988 
Navy tropical global band data. Two filters were applied to 
the data: l) A high-pass filter with a harmonic cutoff at 
20-days, which was applied to remove the strong influence of 
the persistence first guess, and 2) only the first 18 
principal components (EOF's) over a core domain between 
10’S-25"N, 110"-150‘E, were kept for the MCCA. These 18 
EOF's contain approximately 79% of the total variance. The 
two leading MCC modes of the surface meridional wind 
component over the core domain describe a northwestward 
propagating wave pattern with a northeast-southwest tilt, a 
wavelength of approximately 2500-3000 km and a period in the 
range of 6-12 days. These modes appear to be the associated 
with the same wave pattern that was observed by Chang et al 
(1995) in the northern summer 1989-1991 NOGAPS data using 
the same MCCA method. The averaged fractional variance for 
this leading mode pair is about 25% of the first 18 EOF's 
(20% of total high-pass variance), and the average 
correlation between consecutive 12-hour time fields is 0.87. 

Single-year MCCA was performed on each of the 15 
summers and the resulting leading mode structures are 
compared with the 15-year (15Y) modes. In six of the 15 
years(Group A) the single-year mode resembles closely the 
15Y mode with well-defined structure. The composite power 
spectrum of the Group A amplitude time series shows a sharp 
spectral peak at 8 days, which is similar to the 8-9 day 
peak found by Chang et al (1995) . In another five 
years(Group B) the single-year MCC mode resembles the 15Y 
mode partially, with much of the same spatial scale. The 
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Group B composite power spectrum shows a broad range of high 
power between 6-12 days with peaks at 6 and 11 days, 
respectively. The last four years(Group C) have leading MCC 
modes that differ noticeably from the structure of the 15Y 
mode and from each other. The MCC modes computed from each 
of the three groups were compared to the 15Y modes. Group A 
shows the clearest defined wave pattern, Group B also shows 
a pattern that is more clearly defined than the 
individual-year modes of group B years. Even group C modes 
reveal patterns that are similar to the 15Y wave structure. 
Thus, the longer data sets apparently brings out better 
defined structures than the individual year data. 

The wind structure over a larger domain encompassing 
the western Pacific and eastern Indian Ocean between 
20'S-40'N, 7o"E-180 * was then computed using both a 
one-point correlation analysis and a composite analysis 
based on selected cases. Both were done for the 15-year 
data set as well as the three groups. The correlation 
pattern for the 15Y mode at the surface shows that the 
influence of the mode structure covers a larger area than 
the core domain, particularly up and downstream of the 
northwestward propagation track. As was done by Chang et al 
(1995), a search of the locations of tropical cyclone 
centers during times of moderate to large 15Y mode 
amplitudes were conducted and the locations were plotted in 
the correlation maps according to the phases of the leading 
MCC modes. Chang et al. conclud that these tropical cyclone 
centers tend to align along the cyclonic axis of the wave 
mode pattern. The same results were observed for the 15Y, 
and groups A and B modes, and, to a lesser degree, the group 
C modes. 

For the 15Y mode as well as the group A and B modes, 
the wave structure is also reflected, to various degrees, in 
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the 700 and 200 hPa v-correlation maps. The 700 hPa 
structure is nearly in-phase with the surface v, and the 200 
hPa, with lower correlation values, is approximately 1/4- to 
1/2-cycle out of phase with the low levels. For group C the 
correlation at 700 hPa is weaker, and at 200 hPa it is ill 
organized. 

The composite u,v-vectors for 15Y and the three groups 
were done by selecting moderate to strong cases based on the 
amplitude time series, in the same way as the search for 
tropical cyclone locations for the surface v correlation 
plots. The composite was with respect to the time lag from 
the large amplitude time points. For the 15Y and group A 
and B modes, the surface and 700 hPa show couplets of 
cyclonic -anticyclonic circulations that propagate 
northwestward through the core domain. At 200 hPa the 
propagation pattern contains alternating divergence and 
convergence zones that tilt northeast-southwest. These 
patterns indicate that the wave motions are associated with 
strong convective activity, which is consistence with the 
frequent presence of tropical cyclones. Both group A and B 
composites show stronger velocities on the 200 hPa than the 
15Y composite. Group C, as expected, did not show 
consistent organized patterns at 200 hPa. 

The different characteristics between the 11 years of 
groups A and B, and the four years of group C, also 
transpire in the comparison of case selections between the 
15Y modes and the group modes shown in Table 3. The 
majority of the groups A and B cases (62 out of 74) coincide 
(TO within ±12 h) with the 15Y cases (62 out of 70) selected 
in the A-B years. On the other hand, none of the 19 group C 
cases coincide with any of the 20 15Y cases selected in the 
Group C years. It is therefore apparent that the ll A-B 
years possess strong signals of the northwestward 


117 


propagating mode, while the 4 C-years do not. 

A strong correlation between tropical cyclone location 
and the leading 15Y mode was found. Tropical cyclone track 
characteristics are known to vary significantly between 
different flow regimes. The tracks in the northwestern 
Pacific may be classified as straight-movers, recurvers, or 
erratic (Harr and Elsberry, 1995). The climatological 
straight track is more or less along the northwestward 
propagating wave track. It may be conceivable that during 
years of mainly straight tracks, signals of the 
northwestward propagating pattern is stronger, and during 
years of mainly non-straight tracks the signals will be 
weaker. Thus, during years of substantial erratic - track 
tropical cyclone activity the organization of the wave mode 
may be expected to be the weakest. In order to see whether 
the interannual variation in the manifestation of the 
northwestward propagating wave pattern may be related to the 
track characteristics, we compare the tropical cyclone 
forecast errors with the wave mode characteristics each 
year. The hypothesis here is that erratic tropical cyclone 
tracks are the most difficult to forecast. Therefore, large 
72-h forecast error indicates more erratic behavior. 

Table IV lists this comparison where the 72h official 
forecast errors are extracted from the Annual Tropical 
Cyclone Report of JTWC. It is immediately apparent that 
there is no systematic trend for the group A years to have 
smaller errors or group C years to have larger errors. The 
only possible indication of a relation is 1978, which has 
the smallest variance in it single-year leading MCC modes 
(whose 8% variance is 1/3 to 1/5 of the other 14 years' 
values). The 1978's mean 72-h error for typhoons was 459 n. 
mi, which is the largest among the 15 years. The JTWC 
annual report for 1978 also pointed out that it was a year 
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YEAR 

GROUP 

7 2 -HOUR 

ALL TYPHOON 

1974 

B 

348 

357 

1975 

A 

450 

442 

1976 

C 

338 

336 

1977 

A 

407 

390 

1978 

C 

410 

459 

1979 

A 

316 

319 

1980 

B 

389 

362 

1981 

B 

334 

342 

1982 

B 

341 

337 

1983 

B 

405 

384 

1984 

C 

363 

361 

1985 

A 

367 

355 

1986 

A 

394 

403 

1987 

A 

303 

318 

1988 

C 

315 

327 


Table IV:Annual mean forecast errors (NM) 

Western North Pacific 

with unusual, erratic typhoon tracks. However, for 1978 the 
mean 72-h error for all tropical cyclones (410 n. mi) is 
less than 1975 (450 n. mi), which is a group A year. The 
1978 mean error is also comparable to 1977 (407 n. mi), 

which is another group A year. Therefore, the organization 
of the northwestward propagating wave pattern does not seem 
to be related to 72-h mean forecast errors. On the other 
hand, the mean forecast error could be affected 
significantly by a few bad forecasts. In future studies, 
detailed evaluation and comparison of the tropical cyclone 
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track characteristics for the different groups are 
desirable. 

The present data set uses persistence as the first 
guess and therefore its quality is questionable in 
data-sparse areas. New numerical weather prediction 
(NWP)-based, four-dimensional data assimilation datasets 
have the advantage of model forecasts as the first guess, 
but the length of available data sets is relatively short. 
Major NWP centers are now producing reanalyzed data fields 
that start in the late 1970's. Future investigations of the 
northwestward propagating wave pattern may be conducted with 
these improved climatological datasets when they become 
available. In these reanalyzed fields all variables at all 
standard levels are available for a detailed analysis, and 
the results will not be affected by the operational bogus 
procedures. Hopefully, more understanding of the nature of 
the northwestward propagating waves and their relationships 
with the tropical cyclones can be obtained. 
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APPENDIX. 


INDIVIDUAL YEAR MCC MODE #1 STRUCTURE 
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Figure A1: Weighting function for individual years(1974) MCC 
mode#l(equivalent to MCC mode#l structure) of surface 
meridional wind(v) for 12 consecutive 12-hour frames from 
OOh to 132h(5.5 day). Contour interval is 0.06 and dashed 
lines correspond to northerly winds when MCC mode#l 
amplitude is positive. 






Figure A2: Same as Figure A1 except for year 1974 
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Figure A7: Same as Figure A1 except for year 1980 
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Figure A14: Same as Figure A1 except for year 1987 
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